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In this paper we provide an exact formula for the area law of the 't Hooft loop for 
any SU{N) gauge group to cubic order in hot gauge theory. The correction is very 
small for all temperatures above Tc, in stark contrast to the cubic correction to the 
pressure. The gradient approximation in a previous paper, only valid for large N, 
is in excellent agreement with the present, exact evalution. Comparison to lattice 
data is good. Casimir scaling is violated by a small amount not yet resolved by 
the precision of lattice data. 

1 Introduction 

In the study of hot QCD there is a natural quantity, the electric flux loop or 't Hooft loop, 
that measures the electric flux in the plasma. Its behaviour is governed by an area law in 
the hot phase, and, if quark fields are decoupled, a perimeter law in the cold phase. Lattice 
simulations ^ start to determine with considerable accuracy this quantity. 

The coefficient in the area law, the "tension" , is calculable in perturbation theory, in the 
same sense as the pressure ^ . 

Unlike for the pressure the loop is intimately related to a global, discrete symmetry, Z{N) 
symmetry, whose order parameter is the thermal Polyakov loop, winding around the 4d system 
in the periodic temporal or "thermal" direction. 

Like for the pressure there is a loop expansion. But the propagators and vertices in the 
loops arc now propagating in the background of the Polyakov loop. Thus we derive an effective 
potential in terms of the Polyakov loop. 

Z(N) symmetry tells us the effective action is the same for all values of the Polyakov loop 
that differ by a Z(N) phase. In particular the minima of the potential are determined by the 
value of the loop being a Z(N) phase. The minima of the effective potential have the value of 
the pressure. In between the minima the system can tunnel, and this gives rise to the area law 
for the flux loop. 

The pressure has a one loop contribution, the Stefan-Boltzmann result. So has the flux 
loop. The first correction is the two loop result. This is known since long^. The cubic result has 
been derived in the same spirit as for the pressure. The cubic result in the pressure reflects the 
appearance of the Debye mass. In the effective potential the Dcbye mass acquires dependence 
on the Polyakov loop and this gives rise to a natural limitation. The limitation is due to the 
effective potential being concave as function of the Polyakov loop, far enough from the minima. 
This means its second derivative is negative, giving an unstable propagator with a negative 
mass-squared. Nethertheless for a large number of colours one can still establish a result for 
the cubic correction to the flux loop. This was done in the preceding paper ^, henceforth referred 
to as I. This result was in agreement with a recent numerical determination^ of the flux loop 
average in SU(4). 

In this sequel we give explicit results for any number of colours. The method we use is is 
explained in section 2. More details on the method follow in section 3. Predictions are reported 



in section 4, and arc compared with our previous work for large number of colours in the next 
section. The last section concludes. 



2 BVeimework 

The starting point of our investigation is the electric flux loop, as defined in I for SU{N) gauge 
theory: 

VkiL) = eyipi4iT dxdyTrE^Yk. (1) 
Js{L) 

The loop L is taken in the x-y plane, spanning a minimal surface S{L). The electric field 
strength Ek and gauge potential Ak are written in TV x TV matrix form E^ = E^Xa- Quantization 
is obtained through [E^{x), A'^{y)] = \5''^''6k,i6{x - yf . 

The N X N matrix Yk is a generalization of the hypercharge matrix. It generates elements 
of the Z(N) centergroup of SU(N): 

27r 

exp i27rYfe = exp —i—kl (2) 

and can be taken explicitely as: 

Yk = diag(iV -k,N-k, ....TV - k, -k, -k, ... - k). (3) 

To be tracclcss Yk has k entries TV — fc and N — k entries —k. 

Obviously, Vk is an electric flux loop operator and unitary. Because of the canonical 
commutation relations Vk transforms a spatial Wilson loop in the fundamental representation 
by the centergroup factor eq. (2). This property ( the 't Hooft commutation relation^) does 
not depend on the surface S{L). 

Indeed one can show that Vk is a gauge transformation, with a discontinuity expi^ when 
going clockwise around the loop i'*, i.e. the 't Hooft loop^ which is the operator creating a 
magnetic Dirac vortex loop of strength expi^. 

The thermal average of the loop is the Gibbs trace over physical states: 

{Vk{L)) = TrphysVkiL) exp {-^)/Trphys exp (-^). (4) 
In the hot phase the average is an area law: 

(VkiL)) =exp-pkiT)A{L). (5) 

The coefficient Pk(T) is the 't Hooft k-tension. Our interest is mainly in the dependence 
on the strength k. 



2. 1 Pathintegral for the thermal average 

This average can be reexpressed along familiar lines in terms of a path integral. 

As a surface of electric dipoles in the x-y plane our loop will cause, like in static classical 
electrodynamics, a profile for Ao{z) that has a discontinuity across the dipole layer, and dies 
off at infinity. Its slope is continuous everywhere, including at the layer. 

"Normalization of the Gell-Mann matrices is as usual: TrXaXi, = 26"'^ and [Aa, Xt>] = ifabc^c, fabcfa'bc = 

ga,a\ 
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In the non-Abclian case one must find a gauge invariant version of Ao{x). This is the set 
of eigenvalues of the time ordered Polyakov hne: 



P{Ao{x)) = ^PexpiJ^^ dTAo{T,x). (6) 
This is because under a periodic gauge transformation O: 

p{a2) = np{Ao)nl (7) 

Under a gauge transformation Jl; with a discontinuity exp il^ in the periodic time direction 
the eigenvalues will be shifted over Z(N) angles. So in the high T, broken Z(N) symmetry phase 
the Polyakov line is an orderparameter. The eigenvalues will be written as 27rgj, with i = 1, .., N 
and J2i * = 0- 

The average of the loop will then be given by a path integral over all possible profiles: 

{Vk{L))= J DQ{z) j DA(5(exp(i27rQ(^))-P(^o))exp-^5(A) (8) 

where Q is the diagonal NXN matrix with elements g^. The Polyakov line P{Ao) is averaged 

over the transverse directions x and y, as indicated by the bar. 

Doing the path integral over the vector potentials gives the effective action: 



/ 



£>A^(exp (i2TTQ{z)) - P(Ao)) exp = exp-LlSeff{q). (9) 



The effective action has the Z(N) shift symmetries mentioned below eq.(7). 

The transverse size of the system is Ltr- If the loop is very large, comparable to the 
transverse size we can equate A{L) — L^^ in eq. (5). 

As we are interested in the thermodynamic limit Ltr oo the integral over the profile Q 
reduces to minimizing the effective action over the profile with the boundary conditions Q = 
at z = ±oo. The symmetries of the effective action ensure that the minimizing path from 
expj27rQ = 1 to expi^^A; is the straight path from Q = to Q = Yk parametrized by qYk, q 
ranging from to 1. 

This simplifies the task enormously. In terms of the single parameter q we have to minimize 
the efi^ective action with boundary condition q = at z = ±oo and q making a jump of at 
the surface of the loop. 

Finally we find the 't Hooft tension in terms of the effective action by combining eq. (5), 
eq. (9) and the remark underneath the latter to write: 

Pk{T) = minkSeffiq) (10) 

where the r.h.s. is understood to be minimized on the path qYk- 

As explained in the following sections, one can evaluate the tension at high temperature in 
a perturbative expansion, with odd powers of g: 

P^=P^'+pf+pf + .... (11) 

(3) 

The purpose of this paper is to compute p], . 
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2.2 The effective action in perturbation theory 

The effective action can be computed in a loop expansion, in precise analogy with the pertur- 
bative calculation of the pressure. The diagrams are q-dependent through the propagators and 
vertices, due to the delta function constraint in eq. (9). 
Without fluctuations the effective action will just equal 

^^ff(l) = ^^J dzTr{2T,d,qYkf (12) 

in terms of the classical electric field strength E^^ = OzAq. Minimize this with the system 
having extension L in the z-direction. Keeping the kinetic term continuous, whereas q jumps , 
gives a contribution where q{z) = az + b with a = 0(-^). Hence the integral over z will render 

Seff{q)=0{j;). 

So we need quantum fluctuations for the tension. Taking those into account introduces the 
screening length mjj — ^ for (or the Polyakov line). This means that and therefore 
q will vary only appreciably over distances Id = in~j^^. To one loop order the effective action 
becomes the logarithm of the determinant of the fluctuation matrix in the action S{A) in eq. 
(9): 

det{-Dl{q)). (13) 

We suppressed colour indices and gauge parameter dependence, which does not survive the 
determinant anyway. We have Dfj_{q) = 9^ — i5^,o[0/i) and Qq = 2-KTqYk from expanding the 
action S{A) around Q to quadratic order. 

Note that the propagators -^j^ have a background dependence. In the Cartan basis 
propagators with index i ^ j have a mass if (Yk)ii — {yk)jj 7^ 0, because of the commutator 
in the covariant derivative. In all other cases the induced background mass is zero. 

In computing the determinant we can neglect the slow variation of the backgound and get^: 

Seff = -^T j dzTr{2T,d,qYkf + Vi{q) (14) 

with 

Vi{q) = ^^k{N -k)j dzq^{l - \q\)\ (15) 

This potential is periodic in q mod 1 as required by the Z(N) symmetry discussed in the 

previous subsection. 

Upon minimization one gets from eq. (10, (14) and (15): 

as one loop result, with Casimir scaling in the strength of the loop. 
The one loop profile reads: 

, exp^' ,^ , 

qm{z') = 7 17 

1 + exp z' 

with z' = mjjz. 



''Given by the set of N(N-l) generators Xij with only one non-zero element in the (ij) entry, and N-1 
diagonal Ad = (y^2d{d - l))-i/2diag(l, 1, 1, 1 - d, 0, ...0) 
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Figure 1: The only three loop diagram of free energy topology with an infrared divergence. The shaded blob 
is the one loop selfenergy. The colour index of the two propagators need not be the same due to background 
dependence inside the blobs. 



We have, for notational convenience, lifted the branch of the profile on the right of the 
loop by one unit, rendering it continous and ending up at q{z = oo) = 1, not anymore at 
q{z = oo) = . It is for this reason that the tension is often referred to as a domainwall 
tension between two regions of space with different values for the Polyakov line, viz. 1, and 
cxpik^. This is a tenable interpretation of the Euclidean lattice simulation. In a realistic 
plasma however the physical meaning of the Polyakov line is not clear. 

This change in the profile does not alter the minimization procedure because the kinetic 
energy is continuous, and Vi (g) has the required periodicity. The tension does not change. 

The two loop contribution will add g'^V2{q) to the effective action. It will also renormalize 
the kinetic term in the effective action and introduce renormalization effects in the coupling^. 
Quite miraculously Casimir scaling survives in two loops 

2.3 contribution to 't Hooft tension 

Formally the three loop free energy graphs are of 0{g'^). But due to infrared divergencies a 
certain subset of thrc^e loop diagrams is only 0{g^) and this order the only one of interest here. 
They are shown in figure (1). 

In section 5.1 of I we analyzed the order g^ contributions. They originate in the propaga- 
tors sticking into the blobs. If they have no background induced mass then there is a linear 
divergence. 

Note that the same happens for the pressure for all lines. Because the pressure is evaluated 
at q = 0, no line has a background induced mass. All the zero-momentum blobs in this case 
equal the Debye mass rnj^ . And when we integrate the spatial momenta in the graph one finds 
L^{N'^ — 1) J dpj^rn^jj = m^. The point is that all graphs with an arbitrary number of 
selfenergy insertions are of the same order, so have to be resummed. 




Figure 2: Only three loop insertion diagram with infrared divergence 
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This leads for the free energy to a cubic contribution: 

V^ = ^iN'-l)\ogdet{-d^ + ml) (18) 

Returning to the q-dependent case, the value of the selfenergy blob where a line with colour 
a enters, and a line with colour b leaves (both with zero momentum) is denoted by n(A;, q)"''', 
the self-energy matrix. In I we showed this is gauge choice independent. In the limit of g = 
it reduces obviously to m]jSa,b- The indices a and b nm, as said before, through those indices 
that have no background mass induced in the corresponding propagator. In the channel gY^ 
this corresponds to a total number of such indices No{k) = — 1 — 2k{N — k). This is the 
dimension of the self energy matrix. 

The matrix I\.{k,q)a,b can be diagonalized. This is done in I, Appendix C. The result is, 
using the notation r = 

• (k-l)(k+l) eigenvalues equal m|j(l + 6(1 — r)q{q — 1)) 

• (N-k-l)(N-k+l) eigenvalues equal mjj(l + 6rq{q — 1)) 

• one eigenvalue equals ?7i|j(l + 6q{q ~ 1)). 

The total number of eigenvalues corresponds indeed to iVo(fc). 

In the presence of the background q the free energy Vfs changes into: 

Vh + Ft\q) = i(iV2-l)logdet(-52+m|,) (19) 

+ i ( log det{-^5a,b + n(fc, q)a,b) - log det{-^6a,b + n(fc, q = 0)a,b)) ■ 

The last term is subtracted to avoid double counting in the first term. We used a dimensional 
argument that the first term on the r.h.s. of eq. (20), the cubic order of the pressure, is 
proportional to (iV^ — l)L^'m^. 

The second and third term both contain the same No{k)L^m^ contribution, which cancels 
in the difference. This contribution will show up later as the leading contribution in the density 
of states (see the discussion below eq. (24)). But what remains in the difference is 0{L^rnjj) 
i.e. a surface contribution. It is this surface contribution that we will calculate. 

Thus we can put q = qm, because corrections to qm will give higher order corrections to 
the surface term in (20). 

The cubic contribution p^'^-' to the tension is, abbreviating (^m ) = w|j(l+6rgm(<i'm — 1)): 



+ k^ N -k 

+ ^ logdet (- 9^ + m^(g„0) - ^logdet (-9^+ ml,). (20) 

A few comments are in order to transform eq.(20). We work with dimensional regulariza- 
tion, with d — > 3. 

• Use logdet = Tr log and Tr log X - Tr log F = f{exp~tY - exp -tX) 
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• Use the Poisson formula ^„ exp — = -^^i Yli exp — for the periodic trans- 

verse directions. 

A typical term in eq. (20) reads then: 

g{r,d) = -Uogdct{-d^ + ml{q)-hogdct{-d^+ml) 

= iL('^-i) |(47ri)^ (X:exp {-n'LyAt)y'''\Trexp-t{ - + m^iq)) 

-\ Trexp-i(-a2+m|5)|. (21) 

The first trace on the r.h.s. necessitates the evaluation of the eigenvalue spectrum of a one 
dimensional Schroedinger operator: 

-52+m|5(l + 6rg„(g„^-l)). (22) 

3 Solving the Schroedinger equation 

The Schroedinger equation has as potential term 

2 

m^(l + 6rg„((z„ - 1)) = ^(4 - i^rsech'C-)) (23) 
where we used eq. (17). This is a solvable potential^. With 6r = i(i + 1) we have: 

2 

\-dl + ^(4 - j(j + \)sech\\)mz) = E{j),p(z). (24) 

For any integer j we have a bound state entering with parity (— )-' and binding energy 
£'^(j), the superscript referring to parity. There is a continuous spectrum, parametrized by 
E = j)^ mjj. The momentum p runs from to oo. 

The density of states in the positive (negative) parity channel is i/"*" — {v~ — v^). The 
superscript refers to the density of states in the absence of the potential (proportional to ^), 
and should be subtracted according to eq.(21). 

These densities of states can be expressed in terms of the respective phase shifts cos(p.2;-|-^"'") 
{sva.{pz + 5~)). In terms of a dimensionless momentum k defined by p^; = k^^'^ the phase shifts 
are related to the densities by^^^: 

i(l-%]+i,,)(5(fc) (25) 

(26) 

The integer part of j is written as [j], and [j] < j < [j] + 1. Note the delta function in the 
positive parity channel. It is only there when j ^ integer 
The positive parity bound states are for fixed j: 

K^3) = '^{^-{j-nf) (27) 
°We are indebted to R.L. Jaflfe for this reference 



dk 



TT 

°° dk 



{u-{j,k)-VQ) = 



dk.dS+{j, k) 



TT 



TT dk 
dk 85- {j, k) 
TT dk 
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with n a non- negative even integer, and n < [j]. Their negative parity partners are (n odd 
positve integer) 

E^U) = ^{A-{j-nn (28) 

This is important for the correct ultraviolet behaviour of the integral (21) as we will see in 
the next subsection. 

The phase shifts are given explicitely in appendix A. 

We combine the knowledge of the eigenvalues and the relation of density of states to phase 
shifts in a formula for the cubic tension, by plugging eq. (26) into eq. (21): 

g{j,d) = -\{Lmu/^r-'^ l-{i^v)^{Y,e^v{-n\Lrno/Af/v)f-'^ 

X ( cxp -E+ {j)v + Q{j - 1) exp -E^ {j)v 

+ ^°°^(^^-(l-^[,]+i,,)|<5(fc))exp-(fc^ + 4).). (29) 

The 6(a;) = for a; < 0, and 9(x) = 1 for x > 0. The eigentime v = tmj^/A, and the 
bound state energies E are rescaled by the same factor m|,/4 in (27) and (28). 

The phase shift S{j, k) = S^{j, k) + 5~{j, k) and we have absorbed a factor m\,/A into the 
proper time t and called the new variable v. 



3.1 Ultraviolet behaviour 

The ultra-violet behaviour of the tension cannot be worse than that of the pressure in cq. 
(20). And once we have subtracted the ultraviolet divergencies from the temperature induced 
phenomena only infrared infinities may survive. 

The pressure in dimensional regularisation has no ultraviolet infinity. This is best seen by 
writing the pressure in terms of the proper time t: 

1 dt , . -d+l „ 

At t=0 the integrand behaves like and produces a r(— |). Thus dimensional regularization 
subtracts in odd dimensions the u.v. infinities, i.e. produces no poles. 

Let us now look at the expression for a typical contribution to p^. , eq. (29). We get for 
the integrand at v = tmj-,/4: = 0, with Ub the number of bound states: 

(n.-(^^^^ + (l-W.)i)).-. (31) 

This would give for d=3 an ultraviolet divergence. There is, however, a connection between 
the number of bound states and the phase shifts at A; = 0, Levinson's theorem It tells us 
that the coefficient is zero for all j ! 

With the explicit formulae of the previous section we have (j, k = 0) = [j] \ , with 
\j\ = (2n + 1) if j > 2n, n integer. In particular, in between j=0 and j=2, the number of 

positive parity bound states is 1, and ^ ^^'"-^ + ^ = 1. 

For the negative parity phase shifts we have ^ ^''^^ = 9(j — 1) = n^^ for j < 3. 
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More generally the integrand in eq.(29) admits an expansion for small v in powers of w"/^, 
n a positive integer. However for even n the powers in v from the bound states cancel with 
those from the phase shift integral, and only half integral powers stay. This is non-trivial for 
the case of j non-integer, the main reason being that there is not only transmission but also 
reflection. The reader can find the relevant formulae in the appendices. For the convenience of 
the reader we give the formula that results after the cancellation of the integer powers: 

( exp -E+{j)v + e(j - 1) exp -E^ {j)v 

= ^ rduexp-u'smh{{j + l)uV^f^^^^^P^+Io{j,v). (32) 

VTT Jo smh(MVw) 

The remainder Io{j,v) is identically zero, as shown in appendix D. 
The conclusion is that ultraviolet divergencies are absent for any value of j. 
For the special case of j integer there is only reflection. Then, as noticed by Muenster 
one can express the integrand in terms of the error function: 

2 

Erf{^) = -^ dxexp-x^. (33) 
Vi" Jo 

His relation is, for j=l: 

exp {-vE+{j = 1)) + r * ^fciu^ exp -(fc^ + 4)v = exp {-3v)Erf{V^). (34) 
Jo cll^ 

The r.h.s. has an expansion in (v^)^""''"'^) n=0,l,2,... 
Analogously, for j=2 we have : 

exp {-vE+{j = 2)) + exp {-vE^ (j = 2)) + ^^fcM exp -(A:^ + 4)v 

= Erf{2{,/^)+exp{-3v)Erf{y/^). (35) 

The reader can easily check that eq.(34) and (35) are special cases of our general expression 
(32). 

In fact the integral over the proper time v in eq. (29), dropping the L-dependent part in 
the exponent, is elementary in these two cases. For the j=2 case one finds easily, by plugging 
eq. (35) into eq. (29): 

5(r=l,rf = 3)/L2 = _i-(i2 + 31og3). (36) 

For j=l one finds: 

3(r = 1/3, d = 3)7X^ = ^(4 + 3 log 3). (37) 
For interpolating values of j we did the integration numerically (see fig. (3)). 
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3.2 The zero mode at j = 2 
We recall eq. ( 29): 



g{j,d) = f-{nv)^{Y^e^p{-n^Lmn/irMf-'^ 

X ( exp -E+ {j)v + e(j - 1) exp -E^ (j> 

+ f^i^^-^^-^m+^^^m)e^p-{e+4)v). (38) 

As long as the discrete energies are positive (true for < j < 2) we can drop the factor 

( J2n ^^P (~'^^(-^'i^u)^/^)) ^'^ in eQ- (29), because the integral over t is cut-off by the discrete 
energy, so every term in the series is bounded by exp (— n^(Lm£))^£'6/4). That permits us to 
do the large L limit before doing the v integration. 

Hence in the infinite volume limit they are all zero, except the term with n = . And this 
is the way we compute the tension for j <2. 

But for j = 2 we have the positive parity bound state E'^{j = 2) = 4 — = 0. Then any 
term in the series with up to ~ f will contribute, canceling the factor (nv) . For large L 
the result is effectively that obtained in = 1 and we have a logarithmic divergence for v large. 

This is in line with the fact, that in d = 1 wc arc computing the quantum corrections to the 
energy of a classical lump. The uncertainty in the location of the lump at rest is at the origin 
of the divergence in the rest energy. But the zero energy bound state term ("zero mode") is 
needed to give a correct ultra-violet behaviour at u = 0. We can noi just drop the zero mode, 
because then Levinson's theorem would predict an u.v. divergence! 

Muenster introduced a regulator mass, and subtracted in the large volume limit the 
ensuing logarithm. For our purpose we just take as regulator the mass of the positive parity 
bound state. So we do the large L limit before the eigentime integration, and calculate the 
tension as the limit of the integral for j <2. This integral is trivial to do and gives the same 
analytic result as that of Muenster'*. 

For j > 2 the integration over the proper time v starts to diverge because the positive 
parity bound state develops a negative mass. 

4 Predictions 

The basic result is the calculation of g{r,d = 3) in cq.(29). It is shown below in fig. (3). 
Prom the numerical values in this figure and eq. (20), relating g(r,d = 3) to p\. we get 

p'j^'> = (fc2 - 1)3(1 -r,d = 3) + {{N- kf - l)g{r, = 3) + g{l, d = 3). (39) 

In the table below we list for the TV < 8 our prediction for the cubic correction, and our 
old results for lowest order and 0{g^). 

'^The actual value quoted in ref. is diflferent because of finite renormalizations in Muenster's model. In our 
case no finite renormalizations occur (see I). We thank Gemot Muenster for telling us his result in terms of the 

bare variable mo- 
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gHr , d=3L 



0.15 
0.125 
. 1 
0.075 

0.05 
0.025 




Figure 3: The function g{r,d = 3) from eq. (29). 



N 


k 


(a) 


3 


1 


0.370... 


4 


1 


0.598... 




2 


0.769... 


5 


1 


0.831... 




2 


1.190... 


6 


1 


1.066... 




2 


1.624... 




3 


1.797... 


7 


1 


1.301... 




2 


2.068... 




3 


2.422... 


8 


1 


1.538... 




2 


2.519... 




3 


3.063... 




4 


3.236... 



Table 1: Cubic contribution to the tension for AT < 8 
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5 Comparison with gradient expainsion £ind large N results 

Our large N results in I were obtained starting from eq. (20): 

pf = i(A:2-l)||logdet(-5^ + m|,(l + 6(l-r)g(g-l)))-ilogdet(-92+ml,)j 
+ N -k 

+ ^logdet(-a2 + m^(l + 69(9-l))) -ilogdet(-^ + m^). (40) 

We then used the gradient expansion: as in ref. ^ and I we neglect the z-dependence in the 
profile q{z). But then the q-dependent Dcbyc mass in (40) should be non-negative: only those 
terms on the r.h.s. with 1 + 6rq{q — 1) and 1 + 6(1 — r)q{q — 1) non-negative for all values of 
q are allowed. That restricts r = k/N to 1/3 < r < 2/3. These terms are appearing in the 
leading order in a large N expansion, where we keep r fixed, instead of k fixed. 

In that limit the cubic correction to the efi'ective potential becomes: 

= iiv2r2|^logdet(-^ + m^(l + 6(l-r)g(g- 1))) - ilogdet(-a2 + m^)j 

+ r <-^ 1 — r. (41) 

This permits us to do the trace by taking the free particle density of states and doing 
the integral. That gives us a factor i^/(47rt)'^/^. Doing the t-integrations gives us a factor 



m|,r(-3/2) ({I - Qrq{q - l))^/^ - . 
becomes to this order (see ref. I): 



So ' will be of order mjj. The total effective action 



/oo 
dzK(''\q){d,qr + V,^''> + ^ + ^3^=) + ... (42) 
-00 

with V:^''^ as in eq. (49) of I. 

As in I we minimize the action by varying the profile q. The value of the action at the 
minimum q^ is then the tension p^. Doing this leaves us with the result: 

/oo 
dzV^''\qm{mnz)) (43) 
-00 

like in I, eq.(51 and (52)). The only difference is the change of variables from z to q. The 
equation of motion following from the minimization gives us precisely the relation ^ = Vi 
between cq.(43)) above, and (52) in I. 

Note that the minimum profile only depends on the slow variable moz. As the potential 
is 0(m|j) it follows that p^^^ is 0(m|j). 

In fig. (4) we compare the large N result from this paper with the gradient approximation 
in the region where it makes sense. 

As we can see in eq.(39), for large N and finite r (= jj), the leading term in the tension is 

of order A''^. As N goes to infinity, is a function of r, plotted in fig 4. So, in this figure, 

1 ^ 
corrections of 0{j^) have been neglected and only the leading constant term in A'' of is 

plotted. We see that the gradient expansion result is in very good agreement with the exact 

result in the large N limit (There is about half a percent difference.). 
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Figure 4: Gradient expansion (thin line) versus exact computation (bold line) for large JV. 
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Figure 5: ^ up to as a function of . . Dashed line is Casimir scaling, valid for the one loop and two loop 
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prediction. 



6 Conclusions 

The loop can be measured in lattice simulations ^ and in particular ratio's of different strengths 
are relatively easily obtained. Prom the table one can easily infer that Casimir scaling of the one 
and two loop result still holds to a good approximation, even down to the critical temperature. 

Before comparing to lattice data a general comment. Our prediction is in terms of T/Ajjg, 
the lattice data in terms of T/T^. Only for N=2 and 3 the ratio Tc/Ajj^ has been determined. 
So for SU(4) data we take an uncertainty of about 5 % for that ratio. 

In fig. (5) we compare the one loop, two loop and cubic order predictions to the lattice data 
from de Forcrand et al. ^ at various temperatures for the SU(4) ratio P2/ pi- The horizontal 
error bars represent our lack of knowledge of the ratio Tc/Aj^. 

Thus the 't Hooft loop k- ratio is a quantity, that converges well in almost all of the 
deconfined phase for AT > 4. 

There is long standing evidence from lattice data that convergence is good, once the first 
coefficient from the magnetic sector is included. Often this coefficient is the dominant contri- 
bution for reasonable values of T > 2Tc, like for the Debye screening mass^^. As we argued 
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in I, for the 't Hooft loop this contribution will conic into play only at 0{g^), or higher. For 
the ratio's these contributions are apperently not very important. It may be that for the loops 
themselves this correction is needed, to have convergence for T > 2Tc. 
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Appendix A. Phase shifts 

In this appendix the phase shifts in the positive and negative parity channels of the Schroedinger 
equation eq. (22) are given. The reason for this long appendix is that in previous work 
only the solution for r = 1 was discussed. This is an exceptional case where reflection absent. 
In this work we needed all values between < r < 1. Hence reflection is playing a role in the 
discussion. At the same time we have the anomalous relation between density of states 
In terms of C = zm]j/2 eq. (22) becomes 

( - 92 + 6r tanh2(C))V'(C) = (Ar^ + 4)V'(C) (44) 

rrij-) 

and equals q. 12.3.22 of ref. 

It is customary to write 6r = j{j + 1). As j grows, at each integer value of j a bound state 
of parity (— )^ enters. 
Define the shorthand 

e = -^E + 6r - 4. (45) 
rrij-, 

Then, from^ we have for the n'th discrete eigenvalue, n < j: 

en = Qr- (yt + l - (n + i))' = j(j + 1) - (j - n)\ (46) 



It is convenient to use j = '^Qr + j — |. 
Using eq.(45) we get 



m 



En{j) = ^{^-{j-nY). (47) 

We now turn to the continuous spectrum of (73). 

First a few general remarks about parity invariant scattering channels. The wave vector k 
appears in the asymptotic behaviour exp{ikC, + iS^\k)) and exp{—ikQ + i^j[^(fc)) of the phase 
shifts of the transmitted {5^\k)) and the reflected waves {5^R\k)). So the wave vector k is 
related to the momentum p in the phase exp{ipz) hy p = "^"^ . 

We write the wave function with the transmitted and reflected waves as: 

V'^'^ (C) = cxp {ikQ + Rexp{-ikC + (k)) for C « 1 (48) 
^(O(^) = Texp{ikC + id!j.'\k)) ioi C » 1. (49) 

The superscript / indicates the wave travels from the left to the right. R and T are non- 
negative real numbers, the reflection and transmission coefficients respectively. 
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Wc can define a linearly independent solntion with the same wave vector fc, that travels 
from the right to the left. This is just the parity transformation P, with P^' = 

V'^'')(C) = exv{-ikO + Re.xp{ikC, + i5%\k)){ov C,»l (50) 
V;«(C) = Texp{-ikC, + i5P{k)){ov C,«l. (51) 

Alternatively we can define the positive and negative parity continuum states with (we 
drop the r and k dependence in the subsequent formulae): 

^(+)(() = cos(fcC + 5^+^) for C » 1 (52) 
V;(-)(C) = sin(A;C + ^^"^) for C » 1- (53) 

defined to be even (odd) in ^. 

The relation between the two bases is given by: 

+ i exp 

(54) 

and the parity transform of this equation. 

From equations (53) and (49) one finds easily that the reflexion phase shift is related to 
the transmission phase shift by: 

expi((5fl - 5t) = ie(sin(2((5(+) - J^")))) (55) 

and that 

expi{SR + St) = exp {t{d^+^ + S^-'>))te{sm(2{S^+'> - S^-'>))) (56) 

Prom these two equations follows that, as long as sin(2((5(+) — S^~^)) ^ 0, we have the 
identity: 

expi((5(+) + 5^")) = expi^T- (57) 

Only when the positive parity and negative parity phase shifts do lag behind each other by 
multiples of ^ this relation might be invalidated. 

This happens typically when the potential changes from being flat (r = 0) to a shallow 
well (r > O.r << 1. The well then produces only the positive parity bound state. This state is 
formed from a cos{kQ plane wave with k=0. The parity minus sin(fc^) plane waves stay plane 
waves, because the well is not yet deep enough to produce the tanh{Q bound state. But they 
pick up a shorter wavelength due to the well, so are pulled into the well region. 

There is an alternative derivation, simplyfying this condition. 

Prom the definitions of the phase shifts follows: 

T = \T\ expiSx = ^(expi25+ + ex.pi2S-)R = \R\ expiSn = ^{expi25+ - expi2S-) (58) 
Then: 

T±R = cxp2iS±. (59) 
That is: as long as T ^ 0, we factor out T and find: 

2^± =<5T + /mlog(l±|). (60) 
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So 

2{6+ + 6-) = 26T + Im\og{l-{^f) (61) 
and the last term is zero because 

-i^f = -^(m5nsin(2(5+ - 6.)f > 0. (62) 

The case that T = = 1) is only realized for j ^ 1,2 and there only for fc = 0. yields 
for A: = 0: 

± exp iSji = exp 2iS± . (63) 

Hence: 

6+-S.=±^. (64) 
Here the sign depends on the value of j as will be shown in Appendix B. 

7 Appendix B. Density of states in terms of phase shifts 

We put bc's like in Barton 

cos(p„i + 5+) = 0, or p„L + ^+ = 7r(n + 1/2) (65) 
for the positive parity channel and 

sm{pnL + i5_) = 0, or pnL + 5- = nn. (66) 

for the negative parity channel. 

We start with the positive parity channel. We take j ^ 1,2. We find for the sum over 
states with potential, compared to the one without potential: 

e.p-E^{j)t + |^B^g{^-^S{k))e.p{-t'^{k' + 4)). (67) 

For t = we find for this quantity 1 - -1/2 = for j 7^ 1, 2. 

4E+{j)/ml = 4-{jf,j<2. (68) 

For the negative parity channel we have: 

„ , f dkd5-(k), m%,,,o 
e{j-l)exp-E^{j)t + J — ^^(exp-t^(fc2 + 4). (69) 

Here 

4E^{j)/ml = 4 - (i - If, 1 < i < 3. (70) 

For the specific potential we have the phase shifts are now given. For general j we use 
ref.(8): 

T = T{j + l-ik)T{-j -ik)/{T{l-ik)T{-ik)) (71) 

I = r(i-ifc)r(iA:)/(r(j + i)r(-i)) (72) 
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Using Hankels formula, 



r(^)r(i-^) = ^-^, (73) 



to give all arguments in the Gamma funtions a positive real part, one gets for the transmission 
and reflection: 



T{l+j-ik) r(l - ifc)sin(7r(-ifc)) 

sin(7r(-j - iA;))r(l +j + ik) T{l + ik) 



(74) 



ttt/( „:„,_, (75) 



or 



R 

T sin(7rifc) ' ^sin(7r(— j)) 



^ ^(l+j-^fc)^(l + ^fc) sin(7^(-^fc)) 

V{l+j + ik)T{l-ik) sin(7r(-i - ifc)) ' 



(77) 



R sin(7r(— j)) 

T sin(7riA;) 

so for the transmission phase shift: 

(5y ' (fc) = —7m log ( — i sin {Trj)cotanhk + cos(7rj)) 

+ 7m(logr(l +j- ik) - logr(l + j + ik) + logr(l + ik) - logr(l - ik)). {78) 

Binet's formula for the logarithm of the Gamma- functions leaves us with: 



5(:)(fc) = -2 r^Mku) ^\~"7^~'"^^ 

Jo u (exp [u) - 1) 



— Im log ( — z sin njcotanhk + cos nj) 

= B{j,k)+f{j,k) (79) 
where the last equality defines the first and second term. 

As k becomes very large and positive the ratio of the Gamma functions, using Stirlings 
formula in eq. (73) (This is why we transformed to arguments with positive real part.) becomes 
—2Im log(— i) , and the first term Im log exp inj. So the transmission phase vanishes as expected 
for k large and positive. And the same is true for S± from eq.(60). 

The transmission phase is a smooth function of k, apart from a jump in A; = in the term 
/(j, k) = —Im log ( — z sin nj coth nk + cos nj) . 

This function starts for a fixed value of j at large k at the value jn: 

f{h k) = —Im log ( exp —inj + 2i sin(7rj) exp —2k) = nj — Im{2i exp —inj sin(7rj) exp —2k, 

(80) 

so the asymptotic value is corrected by — sin(27rj) exp —2k. 
On the other hand, near k = 0: 

fU:k) = —Imlog[ — isinnj cothnk + cosnj) (81) 

+ -Im\og(-ismnj(^){l + l-{nk)^){l + icot{nj)nk{l-l-{nk)^) (82) 
nk 3 3 

= -Im\og{-i){l + icot{nj)nk) (83) 
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the sign of the coefficient cotg{'Ki)oi k flips at half-integer values of j with period tt. This 
renders the half-integer values oi j fixed points when fc — > 0. 

Hence for < j < 1 the value of /(j, fc = 0) is ^ and 3^ if 1 < j < 2, except in the 
exceptional cases j = integer where there is no dependence of f{j, k) on k and R = Q. It is 
easy to see that: 

/(.7 + l,fc) = /0-,fc)+^ (84) 

so are just rigid translates of one another. This just reflects the growing number of bound 
states and Levinson's theorem, that relates the number of bound states to f{j, k = Q). 

The term ,9± = /m log(l±_R/T) = /mlog (l±i sin(7rj)/ sinh(7rA;)) so equals /mlog(±isin(7rj)) 
at fc = 0. The phase shifts with fixed parity are: 



5±{j, k) = B{j, k) + f{j, k) + g±{j, k) (85) 
Finally, adding up / and g± one finds for < j < 1 : 

<5+(fc = 0) = ^ (86) 

5_{k = Qi) = 0. (87) 

If 1 < i < 2 

<5+(fc = 0) = I (88) 

5_{k = Q) = TT. (89) 

This in accord with Barton's anomalous Levinson theorem for the positive parity channel 
in the presence of only one positive parity bound state, and the normal Levinson theorem for 
the negative parity channel, in the presence of one negative bound state from > 1 on. 

Let us finally look at the exceptional cases, where j is integer. 

Reflection is absent if j = 0, 1, 2 or r = 0, 1/3, 1 so \T\ = 1 there. We have (?± = so 
positve and negative parity phase shifts are identical and equal to the transmission phase shift 
6t- These are precisely the points where a new bound state enters. For r = 1/3, or j = 1, one 
finds from eq.(79), expanding in the argument ku : 

(fc) = TT - 2 arctan k (90) 



and for r = 1 or j = 2: 

5!^^ (k) =2-K — 2 arctan fc — 2 arctan ^ . (91) 

As expected from the general reasoning above, for large k the phase shifts vanish. And 
Levinson's theorem in its normal form is verified: the total number of continuum states, com- 
pared to the number of continuum states without the potential, is 

r '^±5P{k) = 4\k = 0) - S^\k = -oc) ^ ^g^) 

Jo 

with n{r) the number of bound states, from eqns (90) and (91). 



18 



Appendix C. Relation to error function 

Consider eq. (29): 

gij,d) = -i(Lm^/4)(^-i) / ^(7r^;)^(^exp(-n2(Lm^/4)V^))('^-^) 

n 

X ( exp -E+ {j)v + e{j- 1) exp -E^ {j)v 
For j = 1, 2 there is a simple relation to the error function: 



(93) 



2 /"^ 

E{\/v) = / da; exp — a;^ 
Jo 



(94) 



The derivatives of E and of /(?;) = ^prpj exp — + 1) with respect to v equal 
^ } exp —V. That means: 

I{v) = c-E{s/^) (95) 

with c determined by the behaviour at u = 0, c = 1. For j = 1 we have E'^{j = 1) = 3 and 
exp (— 3t;)7(t;), hence 

exp {-vE^{j = 1)) + exp {-3v)I{v) = exp {-3v)E{s/v). (96) 

The r.h.s. has an expansion in (Vw)^"^^, n=0,l,2,... 

For j=2 these same terms are there again, except that oxp— 3t is now furnished by the 
negative parity bound state with energy Ej^{j = 2) = 4 — (2 — 1)^ = 3. We have also a term 
J(t)) = Jp°° ^-ji— exp — t)(A;^ +4), whereas the boundstate with positive parity has become 
zero: 

exp-?;^+(j = 2) = 1. 

Along the same lines: 

J(v) = 1 - E{2^). (97) 

And so for j = 2: 

exp{-vE+{j = 2))+eyip{-vE^{j = 2))+e^p{-3v)I{v) + J{v) = E{2{V^)+e^p{-3v)E{y/^). 

(98) 

Again the r.h.s. is a series in odd powers of ^/v) for small v\ 
For any integer value of j one can easily find the generalization: 

Tr exp -t(9f + V{j, z)) - exp -t(Of + ml) 

j 

= exp-4v(^^eyip{n^v)Erf{Vn?t)y (99) 

n=l 

In this equation the proper time t is related to w by = tm1/4, like before. 
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Appendix D. Interpolation 

For interpolating values of j there is the remarkably simple expression (32) for the trace. It is 
based on the observation that only odd powers of -Jv survive in the error functions for integer 
j and that this generalizes to interpolating j. This we want to prove below. We write the trace 
once more below (for < j < 2): 

( exp -E^{3)v + e(j - 1) exp -El{ji)v 

* - (1 - '5(fc)) exp -(fc^ + 4).) 

2 /""^ , . ,,. ,N ^^ sinh(iMA/u) ^ ,. s , s 

= ^/ duexp-u^sinh i + 1 — ^f^+7oO> • 100 

VTT 7o smh(uVw) 

Remarkable is that the remainder Io(j,v) is identically zero! Use the expressions for the 
transmission in eq. (78) and (79). Integrate the transmission over k. This leads for the 
remainder to: 



smh{uy/v) 



exp {4v)Io{j, v) = exp ((j — n)'^v) j= / du exp —u^cosh{{j + l)u^) — 

- {1-Sjj+,)^- r dk . ,, exp-fc\. (101) 

The last term contains the part of the transmission phase shift f{j,k) defined in (79). It 
vanishes for integer j. So does the one but last term. Obviously only integer powers of v 
appear in the r.h.s. for any j. 

Use now the identity, valid for j integer: 

j 

■ , 7 ^' = exp i-U - l)uV^) V exp 2nuV^. 102 
smh uvt;) 

Substitute this in eq. (101). The second term then reduces to the sum over bound states 
in the first term. So for integer j we have indeed /o(j, v) = 0. 

For non- integer j the coefficients of can be shown to vanish for all I. This proof is related 
to that of the moment sumrules in ref.^^. 

We start by splitting the phase shift B{j, k) in eq. (79) into B{j, k) = Bs{j, k) + Ba{j, k) 
with: 

B,0, t) . r ^isin(,u) ( '^-"■"-^''» + , (103) 

Jo u \^ (1 - exp (-u)) (l-exp(u))y 

and for Ba we have the same but with a minus sign in between the two terms on the r.h.s.. 

After integration over k Bg gives the cosh term in eq. (101) for !(){], v), and Ba the 
sink term in eq.(lOO). We define 5s{j,k) — Bs{j,k) + f{j,k), and Sa{j,k) = Ba{j,k) so that 
6{j, k) = S,{j, k) + 5a{j, k). Note that 5a{j, fc = 0) = 0. 

Then the coefficient of in the expansion of Io{j, v) reads for I > 0: 

[j] 

n=0 
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For I = we get as right hand side in eq. (104): 



1,, . , , Ss{j,k = 0) 
nb-^{l- + , 

which is nothing but Levinson's theorem, eq. (31) so vanishes. 

For / > 1 wc use two facts: 

• -^Ss{j,k) has the poles corresponding to the bound states in the upper half plane (For 
integer j we just proved this.). 

• Ss{j, k) falls of faster than any power of k. 

Thus one can close the contour in the upper half plane, and the coefficient vanishes through 
the cancellation of the sum and the sum of the residues in the integral . 
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